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Soft leptogenesis is a mechanism which generates the matter-antimatter asymmetry of the Uni¬ 
verse via the out-of-equilibrium decays of heavy sneutrinos in which soft supersymmetry breaking 
terms play two important roles: they provide the required CP violation and give rise to the mass 
splitting between otherwise degenerate sneutrino mass eigenstates within a single generation. This 
mechanism is interesting because it can be successful at the lower temperature regime T < 10® GeV 
in which the conflict with the overproduction of gravitinos can possibly be avoided. In earlier works 
the leading CP violation is found to be nonzero only if hnite temperature effects are included. By 
considering generic soft trilinear couplings, we find two interesting consequences: (1) the leading 
CP violation can be nonzero even at zero temperature realizing nonthermal CP violation, and (2) 
the CP violation is sufficient even far away from the resonant regime allowing soft supersymmetry 
breaking parameters to assume natural values at around the TeV scale. We discuss phenomeno¬ 
logical constraints on such scenarios and conclude that the contributions to charged lepton flavor 
violating processes are close to the sensitivities of present and future experiments. 
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I. INTRODUCTION 

Leptogenesis ^ is an attractive mechanism for generating the observed matter-antimatter asymmetry of the Uni¬ 
verse wherein one first creates an asymmetry in the lepton sector which, in turn, induces an asymmetry in the baryon 
sector via anomalous B + L violating interactions. In standard type-I seesaw supersymmetric leptogenesis I2H5] in¬ 
volving the out-of-equilibrium decays of heavy neutrinos and sneutrinos, the CP violation required to generate the 
lepton number asymmetry comes from the neutrino Yukawa couplings. This scenario, with hierarchical right-handed 
neutrinos (RHNs), faces a conflict as successful leptogenesis requires the mass of the lightest RHN to be at least 
10^ GeV m while the simplest resolution of the gravitino problem [710 requires the reheating temperature after 
inflation to be less than 10®“® GeV depending on the gravitino mass (Hln 

One may avoid this conflict by incorporating new elements in leptogenesis. In models of soft leptogenesis nsun] 
(for a recent review, see Ref. |15| 1 CP violation comes from soft supersymmetry (SUSY) breaking terms (here onwards 
we will simply refer to them as soft terms) with soft parameters assumed to be at the tosusy ~ TeV scale; i.e., we still 
hope SUSY is responsible for stabilizing the hierarchy between the weak and grand unification scales. One interesting 
feature is that soft leptogenesis can proceed even with one generation of the RHN chiral superfleld0 Essentially, the 
heavy sneutrino N and antisneutrino N* from the same chiral supermultiplet will mix due to the presence of the 
soft terms. The decays of the mixed mass eigenstates violate both CP and lepton number and generate a matter- 
antimatter asymmetry. Although the CP violation is suppressed by powers of tosusy/AT ^ 1 with M the scale of the 
lightest RHN, the mass splitting between these otherwise degenerate sneutrino mass eigenstates is also proportional 
to tosusy/AT. Grucially, this small splitting also results in enhancement of the CP violation from mixing. Because 
of the suppression factor msusv/A/ in the CP violation, one cannot have very large M. Estimating the leading CP 
parameter as e ^ wsusy/AI and that successful leptogenesis generically requires e > 10“®, we obtain M < 10® GeV 
assuming tosusy Eit the TeV scale. Hence soft leptogenesis occurs in the regime where the conflict with the bound on 
the reheating temperature from gravitino overproduction can be mitigated or even avoided. 

In the original proposals of Refs. mm, the authors showed that in the scenario of V — N* mixing, the leading 
CP violation in decays to fermions and scalars have opposite signs and cancel each other at the order O (uisusy/A/) 
at zero temperature T = 0. They further showed that once finite temperature effects are taken into account, this 
cancellation is partially lifted, i.e. one obtains an asymmetry proportional to a factor \cf{T) — Cb(T)], where cf^b{T) 
are phase space and statistical factors associated with fermion and boson final states, and where the contributions 
do not completely cancel each other at finite temperature. Working under the assumption of proportionality of soft 
trilinear couplings = AY^ where the Ya’s are the neutrino Yukawa couplings and a the lepton flavor index, they 
showed that the resulting CP violation is of the order of O (uisusy/AI) at the resonance which, however, requires an 
unconventionally small soft bilinear coupling B tosusy- Away from the resonance, the CP violation is of O {Y^) 
and, hence, too suppressed for successful leptogenesis. On the other hand, assuming generic A couplings. Ref. [I6] 
showed that successful leptogenesis can be obtained with B ~ tosusy away from the resonant regime. 

Later in Ref. m it was argued that direct CP violation, i.e., from vertex corrections, due to gaugino exchange 
in the loop, survives at the order O (TOlugy/Af^) at T = 0. Since the neutrino Yukawa coupling is replaced by the 
gauge coupling in the CP violation parameter, a large CP violation can be obtained for M at the TeV scale. Further 
study in Ref. m, however, showed that in fact in this scenario, the cancellation still holds up to O (TOg^gy/Af^) at 
T = 0, and it was concluded that finite temperature effects are necessary to prevent the cancellation. The cancellation 
is consistent with the result obtained in Ref. m which states that to have a nonvanishing total CP violation there 
should be lepton number violation to the right of the “cut” in the loop diagram, and this requirement is not fulfilled 
in these cases. More recently, in Ref. it was shown that if finite temperature effects are taken into account 
consistently, the cancellation of direct CP violation from the gaugino contribution still holds even at T 7^ 0. 

In fact, in soft leptogenesis at finite temperature, the partial cancellation in the resulting lepton and slepton number 
density asymmetries sourced hy CP violation from mixing and the complete cancellation in the case of the gaugino 
vertex correction |20| only hold under the assumption of equilibration between the chemical potentials of leptons and 
sleptons (superequilibration) which is valid below T < 10® GeV for tosusy ~ TeV [3. As shown in Ref. m, in the 
nonsuperequilibration regime, the partial cancellation between lepton and slepton number density asymmetries in 
the mixing scenario is avoided, resulting in an enhanced efficiency for soft leptogenesis. However, for reasons given 
later, we shall below consider mixing and vertex scenarios in the superequilibration regime (and also And a case 
where the lepton and slepton number density asymmetries do not partially cancel each other). On the other hand, 
considering M > 10® GeV and tosusy 1 TeV, the CP violating parameter from the gaugino contribution in the 


^ See Refs. I10H12I for another resolution of the gravitino problem due to delayed thermalization of the Universe after inflation. 

^ In a realistic model, we need at least two RHNs to accommodate neutrino oscillations. Assuming RHNs to be hierarchical, soft 
leptogenesis only depends on the parameters related to the lightest RHN and decouples from the parameters related to heavier RHNs. 
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nonsuperequilibration regime is e ~ < 10“^^ and, hence, is too small for successful leptogenesis. 

Therefore processes involving gauginos will not be considered further in this work. 

In this article, we revisit soft leptogenesis by relaxing the assumption of the proportionality of the A couplings. 
In Sec. |TT] we review the Lagrangian for soft leptogenesis with generic Aa terms and spell out the constraints from 
out-of-equilibrium decays of heavy sneutrinos and the cosmological bound on the sum of neutrino masses. In Sec. 

I we obtain the CP violating parameter for both the self-energy corrections (mixing) and vertex corrections. We 
show that generic A^ couplings give rise to two interesting consequences: (1) the leading CP violation can be nonzero 
even when thermal corrections are neglected implying a nonthermal CP violation, and (2) the mixing CP violation 
away from the resonance is of the order of 0{Ya) and, hence, can be large enough for leptogenesis. Because of the 
small mass splitting, the mixing CP violation always dominates over the vertex CP violation even far away from the 
resonant regime. In Ref. [55], it was shown that with A^ = AY^, soft leptogenesis gives negligible contributions to 
the electric dipole moment of charged leptons and charged lepton flavor violating processes. In Sec. |IV[ we repeat the 
exercise and show that with generic A^ couplings, the contributions to charged lepton flavor violating processes are 
close to the sensitivities of present and future experiments. Finally, in Sec. |V| we conclude. This article is completed 
with two appendixes. In Appendix]^ we discuss the inclusion of thermal effects under the assumption of decaying 
heavy sneutrinos at rest. In Appendi x we review the two specific scenarios of A^ discussed in Ref. m and discuss 
an interesting point missed by Ref. |16| which actually allows for nonzero leading CP violation at zero temperature. 


II. THE LAGRANGIAN 

The superpotential for the type-I seesaw is given by 

Wn = + Y,^NfLHu, (1) 

where Nf, £a and Hu denote, respectively, the chiral superflelds of the RHNs, the lepton doublet and the up-type 
Higgs doublet, and i and a are the RHN family and lepton flavor indices, respectively. The SU{2)l contraction 
between £a and Hu is left implicit. In the following, we will assume that the RHNs are hierarchical such that only 
the lightest RHN Ni is relevant for soft leptogenesis. Henceforth, we will drop the family index of RHN, for example, 
N = Ni and Ya = Yiq,. The corresponding soft terms are 

- /:,oft = M^N*N + QsMiViV + Au^m^Hu + H.C.^ . (2) 

The mass and interaction terms involving the sneutrino N from Wn are given by 

-Cj^ = |M|^ N*N + (^M*Y^N*I^Hu + Y^H^Pl£o.N + H.c.) , (3) 

where Pl,r = | (1 T Ts)- Through held redefinitions, it can be shown that the three physical phases are 

$„ = arg(A„F„*R*). (4) 

Without loss of generality, the phases can be assigned to Aa and all other parameters will be taken real and positive. 
We would like to stress that we do not assume the proportionality of Aa to the neutrino Yukawa couplings (Aa = AYa) 
as has been done in Refs. |I3||I11|53| where there is only one physical phase $ = arg(AR*). As we will show in Sec. 
|HI| by considering generic Aa couplings, the CP violation can be nonvanishing even at zero temperature. 

Because of the bilinear B term, N and N* mix to form mass eigenstates 

n+ = ^(n + n*), 

N_=-^(n-N*), (5) 

with the corresponding masses given by 

Ml = + M^± BM. (6) 

In order to avoid a tachyonic mass which implies an instability of the vacuum such that the sneutrino will develop a 
vacuum expectation value, we always assume B < M + 
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Rewriting the Lagrangian in terms of mass eigenstates N± we have 


-Cj^- Aoft = MlN^N+ + MlN*_N_ 




+ (A^ + MY^) 


+iN_ \Y^H^PL^c. + (^a - + H.c.} . 


(7) 


A. General constraints 


The total decay width for N± is given by 


r± 


M y2 \Aa\'^ yQRe(AQ,) 
“ + 2M2 M 

a *- 


( 8 ) 


where we have expanded up to 0{Y^, /M'^,Yams\jSY/AI) and ignored the final state phase space factors. We 

will impose the restriction that \Aa\,B < M and < 1 to ensure that we are always in the perturbative regime. In 
principle, msvSY/M and Ya can go up to 47r before perturbative theory breaks down but with our stronger restriction, 
we are not anywhere near the nonperturbative regime. 

The out-of-equilibrium condition for leptogenesis is 


r± <H{T = M), 


(9) 


where the Hubble expansion rate is given by 77 = 1.66y^ /Mp\ with Planck mass Mpi = 1.22 x 10^® GeV. Assuming 
minimal supersymmetric Standard Model relativistic degrees of freedom, we have 5* = 228.75. The condition above 
translates to 


Y2 + 


|A 


_" YJie[A^) 

2M2 M 


< 1.6 X 10 


-5 


/ M 

107 GeV ) 


( 10 ) 


From the condition above, we see that \Aa\ is bounded from above depending on M. For example if M ^ TeV, we 
require \Aa\ < 10“^ GeV. At this low scale, the mass splitting between A^_|_ and iV_ is required to be of the order 
of their decay widths such that the CP violation is resonantly enhanced to yield successful leptogenesis |25l l26] . To 
avoid excessive fine-tuning, if we consider \Aa \ ~ TeV, Eq. (101 implies M > Ax 10^ GeV. 

In type-I seesaw, barring special cancellation, we have the upper bound on the sum of light neutrino masses from 
cosmology fM\ 

w2 2 

~ 0.23 eV, 

M ^ ' 




V? < 3 X 10 


-4 


M 


107 GeV 


1/2 


1-f 


tan^ j3 


1/2 


( 11 ) 


where tan,d = Vujvd and Vu{d) = {Hu{d)) are the up(down)-type Higgs vacuum expectation values, -I- wj = 
y/2 ~ (174 GeV)^ with Gp the Fermi constant. For tan ,5 > 1, the bound above is always less stringent than 

Eq. (|10[), and, hence, the out-of-equilibrium condition alone is sufficient. 


III. CP VIOLATION 


In this section we will study CP violation of the Lagrangian ([^ from the interferences between tree-level and one- 
loop diagrams shown in Figs. [l]and[^ We will take into account thermal corrections while approximating sneutrinos 
N± to always be at rest with respect to the thermal bath. Since we are in the regime where all three lepton flavors 
can be distinguished (T < 10® GeV), we will not sum over the lepton flavor in the final states |23| . 

To quantify the CP violation, we define the CP asymmetry for the decays N± —>■ Oq, with Oq, = {£aHu,iaHu} as 


S,v ^ 

^iba — 


7(iV± o„) - 7(7V± a„) 


Ea, 3;/3 7(^± ^ a^) + 7(7V± ^ a/d) 


( 12 ) 
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Figure 1. One-loop self-energy diagrams for the decays N± —^ iaHu [(a),(d)] and N± —> £aHu [(b),(c),(e),(f)]. The arrow 
indicates the flow of lepton number. The dotted vertical lines indicate the corresponding intermediate states go on mass shell. 
The diagram with fermionic loop and fermionic final states does not contribute to the CP violation since it does not involve 
the soft couplings Ac- 


where the superscripts S and V indicate the CP violation coming from self-energy and vertex corrections, respectively, 
cCi indicates the CP conjugate of Ua, and 7(1 —>■ j) is the thermal averaged reaction density for the process i ^ j 
defined in Eq. ( |A1[ ). In the following, we will include the thermal effects associated with intermediate on-shell states 
which, as shown in Ref. |20| . will result in the cancellation of vertex CP asymmetries from gaugino contributions 
nziiiH]. We will always approximate N± to be at rest with respect to the thermal bath so that we can obtain analytical 
expressions for the CP asymmetries (see Appendix]^. Furthermore, we focus on the superequilibration regime which 
falls in the temperature range T < 10® GeV for msusY ~ 1 TeV j5j. The advantage is that in this regime, lepton and 
sleptons are not distinguished (they have the same chemical potentials) and so the two Boltzmann equations for the 
lepton asymmetry in particles and sparticles can be reduced to one equation for the net lepton asymmetry]^ Hence 
we are allowed to sum over CP asymmetries of lepton and slepton final states as below. 


A. CP violation from mixing 


In this subsection, we discuss the mixing CP violation from self-energy corrections. There are two kinds of self¬ 
energy diagrams as shown in Fig. the diagrams with continuous flow of lepton number [Figs. l(a)--l(c)] and the 
diagrams with flow of lepton number inverted in the loop [Figs. l(d)-l(f)[. Notice that diagrams with fermionic loop 
and fermionic final states do not contribute to the CP violation since they do not involve the soft couplings Aa- From 
Figs. l(a)-l(c), we obtain the respective contributions to the CP asymmetries defined in Eq. (12) as followffl 


S,(a 

e±a 


S,(b) 

e±a 


S,(c) 

£±q 


1 

47rG±(r) 


Im(A/3) 


B 


2BM 


AB^ 


rB{T)cF{T), 



(13) 


2BM 
' 4R2 + r| 


rB(T)cB(T), 


^ We make the assumption of superequilibration also to highlight the positive effects of nonthermal CP violation in soft leptogenesis. 
Including nonsuperequilibration effects, the efficiency of soft leptogensis is expected to be further enhanced, and this effect was studied 
in detail in Ref. [ 21 ]. The validity window of superequilibration can be enlarged by increasing the gaugino masses and /r parameter [5] 
and/or decreasing |Aq:|. 

^ The absorptive parts which regularize the singularity in the N± propagators as M+ M— are obtained by resumming self-energy 
diagrams following Refs. |25II26| . 
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where we define 

G±{T) = 


y2 _|_ ^ ^ 2yaRe(^a) 


M2 


M 


c.(T) + y^ l + ^±^ c.(T). 


(14) 


In the above TB^piT) and Cb.f{T) are temperature-dependent terms associated with intermediate on-shell and final 
states respectively, as given in Appendix We will also make use of the following identity 


rF{T)cB{T)=rBiT)cFiT), 


(15) 


proven in Appendix 


Note that if we sum over the lepton flavor a and use Eq. (151, we obtain + £±1^^) = 

Ea = 0, in agreement with the T — 0 result of Ref. [H] that if there is no L violation to the right of the cut in 
the one-loop diagrams, the net CP violation on summing over all final states is zero. 

From Figs. l(d)-l(f), we have 


S,(d) _ 

^±q; 


47rG±(T) 




Ini(Aff) ‘iBM 

M 1 M2 M ) 4R2 


S,(e) _ 

e±a — 


SXf) 

^±a 


AfG±{T) 

1 

47rG±(r) 


V2v 

“ M 


, r\ 2BM 

^ M^ M ) 4R2 +r2 


(16) 


-k^-E 


\ Im(A.) / |A 
M I “ M2 




Im(A/3) 


M 


2BM , , , , 

^ 4^2+r2 


Notice the leading contributions from Ab|_ and 7V_ in Eqs. (131 and (161 come with the same sign and, hence, they 


will contribute constructively to the lepton number asymmetry. 

The total CP asymmetry from mixing = J2n={a,b,c,d,ej} given by 


4„ = 


Im(A,g) ABM 


^ y2 Y 

AnCiiT) ^ 

'^47rG±(r) M'^ ^ M 4R2+rE^ 

1 Im(A«) / M2 b\ ABM 


(17) 


47rG±(r) 


“ V M I M2 M j AB^+Tl 


rB{T)cF{T). 


In the above, the first term vanishes in the zero temperature limit T —>■ 0 when cb,f{T) —> 1 and TB^piT) —>■ 1, while 
the terms higher order in msusv/M survive. They remain nonzero after summing over the lepton flavor a. In the 
following in order to make the dependence of thermal and nonthermal CP asymmetries in Eq. ( [T7| on the model 
parameters more transparent, it is instructive to look at two limiting cases (i) Ya ^ A^/M and (u) Y^ ^ Aa/M 
where in case (i), the thermal CP violation dominates, while in case (ii), the nonthermal CP violation dominates 


In the limit (i) Y^ A^/M, we have 

1 


4o 


__Lp Im(2^/3) 4RM cf{T)-cb{T) 

~ Af^ /5 ^ AB‘^+TI.cf{T) + cb{T) 


(18) 


where we define the flavor projector P^ = Y^jY'^ with Pa = ^ and T^ = , and we have dropped the 


terms higher order in msusv/Af. In this case, the CP asymmetry in Eq. (18) is proportional to cf{T) — cb(T) 
which goes to zero as T —0, and, hence, the contribution to the CP violation is the thermal one. 


® By thermal (nonthermal) CP violation, we refer to the case where CP violation does (not) vanish as T —> 0. 
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In the resonant regime where B ~ r±, we have e|. ~ {\A\/M)/Y where we have suppressed the lepton flavor 
index for an order of magnitude estimation. In this case, a large ef. can be obtained which allows TeV-scale 
leptogenesis but at the cost of having unnaturally small ^ TeV. 


Away from the resonant regime when B ^ r±, the CP asymmetries go as ef ~ 10“^y|A|/i? assuming 0{1) 
contribution from the CP phases of Eq. Q. Taking |A| ^ TeV > B together with the out-of-equilibrium decay 
condition (101 gives us sufflcient CP asymmetries e|. > 10“® for M > 10^ GeV. 


• In the other limit (ii) Ya Aa/M, we have 




1 4i3M 


(19) 


fvi! 

a S-ttM 


The CP asymmetries Eq. (191 clearly do not vanish at T = 0, and this represents 


where Tyi = ^ 

a nonthermal CP violation. Of course thermal effects are always there but the fact that the CP violation is 
nonvanishing at T = 0 implies that it is less suppressed compared to case (i). 

In the resonant regime B ~ r±, we have ~ Y/{\A\/M). In this case too a large can be obtained which 
allows TeV-scale leptogenesis but at the cost of having unnaturally small |A|,i3 <C TeV. 

Away from the resonant regime with B ^ r±, the CP asymmetries, like in the limit (i), go as ~ 10“^V|A|/i? 
assuming 0(1) contribution from the CP phases of Eq. Q. Hence taking |A| ~ TeV > B together with the 
out-of-equilibrium decay condition (10) gives us sufficient CP asymmetries > 10“® for M > 10^ GeV. 


To confirm our estimation of successful leptogenesis and also to illustrate the enhancing effects of nonthermal CP 
violation, we numerically solve the Boltzmann equations using the expression for the asymmetry parameter in Eq. (171. 
For simplicity, we consider only decays and inverse decays of N and N±. We will also define the washout parameter 
as K = r±/H{T = M) with r± given by Eq. (| 8 1 . In Fig. we plot the absolute value of the final baryon asymmetry 
|Vab(oo)| as a function of K for the following three scenarios: 

• Non Thermal dominated (NTD): In this scenario, we choose A/M = (10“"^, 10“^, l)ri; and V= (10“®, 10“^, 10“^)ri 

• Thermal Dominated (TD): In this scenario, we choose A/M = (10“®, 10“^, 10“^)w and Y= (10“^, 10“^, l)w. 


• Mixed (MIX): In this scenario, we choose A/M = (10 10 l)w and Y= (10 10 l)w. 


In the above A and Y represent the couplings written as 3-vectors. In all the scenarios above, we vary w between 10“® 
and 10“"^ such that we scan through the parameter space from the weak washout {K = 0.1) to the strong washout 
{K = 15) regime while still respecting the cosmological bound on the sum of light neutrino masses in Eq. (111. For 
definiteness, we also fix M = 5 x 10^ GeV, tan/3 = 10, arg(AQ,) = —7r/2, and B = 1 TeV. 

In Fig. 1 ^ for the left plot, we solve from an initial time with T M assuming zero initial number densities for 


N and N± while for the right plot, we assume thermal initial number densities for N and N±. For the observed 
baryon asymmetry, we use the recent combined Planck and WMAP GMB measurements of cosmic baryon asymmetry 
[M HZ] at 2a, 


yCMB ^ (-g gg ^ Q 22) X 10"“ , 


( 20 ) 


which is plotted as the gray band in Fig. From the plots, we see that in the NTD (blue dashed line) and MIX 
(purple solid line) scenarios, |yAB(oo)| falls off very slowly in the strong washout regime AT > 1. The reason is that 
the falloff in their efficiencies is almost completely compensated by the increases in their respective CP asymmetries 
as one increases Aa/M. On the other hand, in the TD (red dotted line) scenario, in the AT > 1 regime, |Vab(oo)| falls 
off much faster due to the additional suppression from the partial cancellation between the CP asymmetries from the 
decays of N± to scalars and fermions. 

In our study we are also interested in the situation when Aa and B are restricted to be around the TeV scale. In 
Fig.H the regions to the left of the thick blue dashed and purple solid vertical lines correspond to AT when Aq, < 5 
TeV for the NTD and MIX scenarios, respectively, while Aa < 5 TeV for the TD scenario in the entire range of AT 
considered in the plot. We find, for example, for the case of zero initial number densities of /V and JV±, the correct 
amount of baryon asymmetry can be obtained for NTD at AT ~ 0.8, MIX at AT ~ 0.6, and TD at AT ~ 4. Notice that 
the appropriate sign baryon asymmetry can always be obtained by choosing the appropriate phases of the complex 
couplings Aa- From this numerical exercise we conclude that generation of sufficient baryon asymmetry is possible 
for TeV-scale Aa and 73 ^ r±, i.e., far away from the resonant regime. Besides, we also see that nonthermal CP 
violation can significantly enhance the efficiency of soft leptogenesis. 

Finally, in Appendix^we will discuss two special cases, namely, (a) Aa = AYa and (b) A„ = AY"^/(iYa) considered 
in previous work. 











5 = 1 TeV, arg(Att) = -:7r/2, tanjS=10 


5 = 1 TeV, arg(A(j) = -:7r/2, 


tanjS=10 


IFasCoo)! 






Figure 2. The absolute value of the final baryon asymmetry |Yas(oo)| as a function of the washout parameter K = V±/H{T = 
M) for M = 5 X 10^ GeV for the three scenarios described in the text: NTD (blue dashed), TD (red dotted) and MIX (purple 
solid). The left plot corresponds to the case of zero initial number densities of N and N±, while the right plot corresponds to 
the case of thermal initial nnmber densities of N and N±. The regions to the left of the blue dashed and purple solid vertical 
lines correspond to K values when Aq, < 5 TeV for the NTD and MIX scenarios, respectively, while for the TD scenario we 
always have Aq, < 5 TeV in the range of the plot. The gray band represents the recent combined Planck and WMAP CMB 
measurements of cosmic baryon asymmetry |241127| at 2cr. The dip in the TD scenario in the left plot refers to a change in the 
sign of the baryon asymmetry. 



Figure 3. One-loop vertex diagrams for the decays N± —>■ laHu [(a)] and N± —>■ icHu [(b),(c)] with the conventions of Fig. 

The diagrams with fermionic loop and fermionic final states do not contribute to the CP violation since they do not involve 
the soft couplings Aq. 


B. CP violation from vertex corrections 


In this subsection, we discuss the CP violation from vertex corrections. From Figs. 3(a)-3(c), we obtain 


V,{a) 

e±a 


V,ib) 

^±a 


V,{c) 

e±a 


T 


T 


± 


1 


87rG±(T) 


M - M~2 - rB{T)cF{T), 


M 


1 


87rG±(T) 

1 




Im(AQ) Ml + M2 


M 


87rG±(T) 

m2 Mi + Ml 

X —^ In-^ 

Ml Ml 






M 2 


rF{T)cB{T), 


M 2 
P 

^tb{T)cb{T). 


Yr 


Im(AQ,) 

M 


+ Y^- 


1^0 


M2 




Im(A/ 3 ) 

M 


( 21 ) 
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Figure 4. One-loop vertex diagrams for the decays N —>■ icHu [(a)] and N —>■ icHu [(b)] with the conventions 


of Fig. 


Summing over the contributions above and expanding in B/M ^ 1 in the numerators, we have 

(c) 


_ V,(a) V,ib) v,(c) 
^±a = e±a + £±a + ^±a 


= T 


In 2 

87rG±(T) 

1 

'87rG±(T) 

In 2 

■ 87rG±(T) 

1 

' 87rG±(T) 




M 


+ YaY.^P 


M 


[cpiT) - CB{T)]rB{T) 


Im(Aa) 
“ M 


+ 

' ^ Oi 




B 

M 


cf{T) 


(In2-1)CB(T) 


rsiT) 


( 22 ) 


E 

. (9 

E 

(9 


\A^ Im(Aa) 

M 2 “ M 


\A^ Im(Ac,) 

M 2 “ M 


M2 

M2 






M 

Im(A/ 3 ) 

M 


rB(T)cB(T) 


M 


(ln2-l)rB(r)cB(r). 


The leading contributions from 7V± [first and third lines of Eq. (22l] come at the order of ^ 10 [taking 
Yq, ~ Im(AQ,)/M], which are too small for successful leptogenesis from Eq. (101 for M ^ 10^ GeV. Of course the same 


conclusion holds also when Yq, ^ lm(^Q,)/M or Yq, <C lm(AQ)/M. Besides, notice also that the leading contributions 
from N± come with the opposite signs, and, hence, they will contribute destructively to the total lepton number 
asymmetry. Upon expanding G±{T) terms also in B/M ^ 1, we obtain an additional suppression factor B/M like 
the terms in the second and fourth lines in Eq. ( |2^ . Hence we conclude that the vertex CP violation is irrelevant for 
soft leptogenesis. 

So far, we have been discussing the contributions of soft terms to CP violation in the decay of N±. In fact the soft 
terms also provide new sources of GP violation in the one-loop vertex diagrams for the decays of the heavy neutrino 
N as shown in Fig. Nevertheless the CP violation from these diagrams comes at the same order as Eq. (221 and, 
hence, is too small for successful leptogenesis. 


IV. PHENOMENOLOGICAL CONSTRAINTS 

We are primarily concerned with scenarios with M 10^ GeV for which the production of sneutrinos is beyond 
the energy range of current colliders. However, even if M± ~ TeV, the bound on the Yukawa couplings from the 
requirement of out-of-equilibrium decays of N± [in Eq. (H^i makes N± impossible to be produced at colliders |28| . 
On the other hand, the soft SUSY breaking parameters relevant for soft leptogenesis Aa, B, and M can contribute 
to Electric Dipole Moments (EDM) of leptons and to Gharged Lepton Flavor Violating (GLFV) interactions though 
the analysis of Ref. [22] under the assumption of universality soft trilinear couplings Aa = AYa showed that the 
contributions to EDM and GLFV are much below the experimental bounds. Here we will repeat the analysis of 
Ref. [55] considering a generic Aa- In Ref. [IS], the phenomenological consequences of the soft terms considering three 
generations of RHN chiral superfields have been discussed at length. Glearly, these soft parameters are connected with 
the mechanism of SUSY breaking and as such are model dependent. Here we will remain agnostic about the SUSY 
breaking mechanism and simply focus on the phenomenological constraints on these parameters and, in particular, 
we will focus only on parameters related to 7Vi which are relevant for soft leptogenesis, i.e., B, M, Aq, Yq,, and M. 
Without fine-tuning, we consider the soft parameters B, M, and Aa to be similar or smaller than tosusy TeV. On 
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the other hand, the parameters Yq, and M are subject only to the out-of-equilibrium N± decay constraint in Eq. (101 
and less stringently to the cosmological bound on the sum of neutrino masses in Eq. (11). The running of Yq from 
the high scale down to the weak scale gives some corrections at the level of 10% — 20% (see Fig. 3 of Ref. [3T]) which 
we will ignore in the following. 


1. Electric dipole moment of the electron 


Assuming 0(1) contribution of the phases and mixing angles in the chargino sectors, the contributions of Aa and 
B to the EDM of the electron are given by 


\de 


e rrif, tan ^ 
IOttto? 


(14 

M2 ° 


BY^) , 


(23) 


where me is the electron mass, is the squared mass of the light sneutrino and the mass of chargino. For generic 
Aa, the first term in in Eq. (23) dominates. Taking mp = = m-suSY and making use of Eq. (10), we have 

^_ 3 o /tan/3\ /lO^GeVV^^ f ITeV \ 


\de\ < 5 X 10 


V 10 


M 


■ 


\msijSY J 


e cm, 


(24) 


which is much stronger than the current experimental bound |de|exp < 8.7 x 10“^®ecm [30] . The contributions to /i 
and T EDM can be estimated by replacing mg in Eq. (23) by and m^, respectively, but the current experimental 
constraints on them are still a lot weaker: |(ip|exp < 1-9 x 10“^®ecm |32] and |c?r|exp < 5.1 x 10“^^ecm |33| . 


2. Charged lepton flavor violating interactions 

The branching ratio for charged lepton flavor violations due to nonvanishing off-diagonal elements of the soft mass 
matrix of the doublet sleptons to? is given by |22l [34] 


BR{ia ^ ihl) 


a 


(mj) 


"^iuSY 


tan^ /3, 


(25) 


where a is the fine structure constant. In general, the off-diagonal elements of to| will induce too-large CLFV rates. 
The usual solution is to assume mSUGRA boundary conditions at the grand unified theory (GUT) scale where the 
off-diagonal elements of to? vanish. In this case, as to? evolves from the GUT scale Mqut to the RHN mass scale M, 
the off-diagonal elements will be generated due to the renormalization effects as |3S] 


(w^)o!/3 




/ Mgut \ 

) 


(26) 


for a 7 ? /?, and we have kept only the dominant contributions from Aa- 

The most stringent constraint on the rare decay /i —>■ ey comes from the nonobservation of the process from the 
MEG experiment IM1I371 which has set the new bound on the branching ratio for fi —> ey. 


BR(/r —>• ey)exp < 5.7 x 10 


Substituting Eq. (26) in Eq. (25) and applying the constraint Eq. (27), we obtain 


< 5 X 10^ GeV^ 


/ WSUSY 

V ITeV 



(27) 


(28) 


where we have taken Mgut = 10^® GeV and M = 10^ GeV. Similarly using the experimental bounds on GLFV in r 
decays, BR(r —>• ey)exp < 3.3 x 10“® and BR(r —)• ^y)exp < 4.4 x 10“® j3H]> we obtain 




|A*A^| < 1 X lO^GeV^ 


/ nT-SUSY 
V ITeV 



(29) 
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For TOsusY at the TeV scale, the bound (29) can be satisfied with Aa also at the TeV scale while the stronger bound 


(281 requires either A^, and/or A^ to be smaller than TeV scale. As discussed in Sec. Ill A the mixing CP asymmetries 
away from the resonant regime go as ef ^ 10“^yQ,|AQ,|/i? and can be large enoug^h with M > 10^ GeV and having 
one of the A^ ~ TeV > B. 

In addition, the off-diagonal entries of the slepton mass matrix can also give rise to other CLFV processes like 
/i —>■ 3e and /i — e conversion. If such processes are dominated by the dipole-type operator for relatively large tan /3, 
BR(/i —> 3e) and the rate of — e conversion rate R^e ^-re proportional to BR(/i —> ey) and are approximately given 
by m 


and, for the fjjAl nucleus, by |1D] 


BR(^ —>• 3e) ^ 6.6 x 10 '^BR(^ —>■ ey) 


Rfie ~ 2.5 X 10 '^BR(/i — ey). 


(30) 

(31) 


The present constraints coming from these CLFV processes are less severe than those coming from p. —> ey. However, 
in future experiments the sensitivity for such processes may improve, which could constrain the presently allowed 
parameter space or lead to a detection of such lepton flavor violating processes. As for example, the future Mu3e 

experiment im could reach a sensitivity of ^10 -10 for BR(/i —>■ 3e). For the /r — e conversion process, from 

the Mu2e |32] and COMET |33] experiments, the bound could reach the level of ~ 10“^’’ for the f^Al nucleus, 
while the PRISM/PRIME [IS] project may have 2 orders of magnitude greater sensitivity. 


V. CONCLUSIONS 

In the framework of local SUSY, soft leptogenesis is an attractive mechanism to explain the cosmological matter- 
antimatter asymmetry since it works at the lower temperature regime T < 10® GeV where the conflict with the 
overproduction of gravitinos can be relaxed or even evaded. We showed that by considering generic soft trilinear 
Aa couplings there are two interesting consequences: (1) one can realize nonthermal CP violation where the CP 
asymmetries in the decays of heavy sneutrinos to lepton and sleptons do not cancel at zero temperature resulting in an 
enhanced efficiency in generating baryon asymmetry, and (2) the dominant CP violation from self-energy corrections 
is sufficient even far away from the resonant regime and the relevant soft parameters can assume natural values at 
around the TeV scale. For successful soft leptogenesis, we considered two requirements: the out-of-equilibrium decays 
of heavy sneutrinos and a large enough CP violation. Assuming tosusy TeV, we found the following conditions 
Aa ~ TeV > B and M > 10^ GeV as sufficient for successful leptogenesis. In addition we also found that while the 
contributions to the EDM of charged leptons are negligibly small, the contributions to CLFV processes are close to 
the sensitivities of present and future experiments. 
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Appendix A: Thermal corrections 

The thermal averaged reaction density is defined as 


y {ab... ij...) = A))/;; \M {ab... ^ /^/r- (1 + VifD (l + ■■■. 


(Al) 


where Ai{ab... —t ij...) is the amplitude for the process ab... ij... at finite temperature, = (e^’/^ — rji) ^ with 

rji = ± for i representing the boson or fermion, respectively, and 


Ktj... = 


j (niadflb...dflidl\.j... 5^"^^ {pa +Pb + ... -Pi -Pj - •■•), 


dHj = 


(Ppi 

( 2 ^)" 2E, ■ 


(A2) 
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The CP asymmetry for the decay a —>■ ij which arises from the interferences between tree-level and one-loop 
diagrams as shown in Fig. [^is defined as 


7 (q -» ij) -j{a^ ij) 
T.k,i [l {a ^ kl) +-f {a ^ kl)] 


f (1 + (1 + rj.fp) 

\M (a ijp - 

M [a ^ ij)\^ 


Ek,l I dUafa^ f (1 + VkfP) (1 + Vi/n 

\M{a^kl)f 

+ 

M(a^ ki) 

|2' 


where the two-body phase space integral is 

J d^ij = j dUidlij {2Tr f {pa - Pi - Pj) ■ (A4) 


Ignoring the thermal motion of a with respect to the thermal bath, i.e. setting Ea = Ma, we can drop the integral 
JdU, in both the numerator and denominator. In this case the phase space integral can be carried out analytically. In 
order to obtain the thermal factors associated with the intermediate on-shell states, one necessarily needs to calculate 
the amplitudes in Eq. (A3) using thermal field theory. In the real-time formalism of thermal field theory, one needs to 
double the number of degrees of freedom (introducing type-1 and type-2 fields) resulting in a 2 x 2 matrix structure for 
the thermal propagator (see, e.g., |SI])- However, at one loop, we can take all the vertices connected to external legs 
to be of type 1 and, hence, we only need to consider the (11) element of the thermal propagator. The (11) component 
for the boson propagator is m 


^- TnTT- + - msirn, (A5) 

p^ — niByd ) + 

where msiT) is the boson thermal mass, and the cut propagators are 

D% = 2tt [e{±po) + /sdpoD] 5{p'^ - mB{Tf). (A6) 

For fermions, the structure of the propagator is more involved m- For simplicity, we approximate the (11) fermion 
propagator by 




- 27r/i^(|po|)^(p^ - mpiTf) 


_ rnp{TY -I- ie 

where mpiT) is the fermion thermal mass and the cut propagators are 

D$ = ^TTf [e{±po) - fpilPol)] ^(/ - mpiTf). 


(A7) 


(AS) 


In Eqs. (A7 \ and (A 81 , the propagators ~ p are without a mass term as the bare fermion mass is zero, and the thermal 


mass does not have chiral properties. Also, as implicit in the propagators, we have considered the dispersion relation 
as that of a free particle with a thermal mass, instead of the actual dispersion relation including thermal corrections. 
Although this is an underestimate of the actual dispersion relation, the error is within 10% jH]. The above also 
implies that in Eq. (A7l we have ignored the fact that due to the interactions with the thermal bath the two poles 


of the fermion propagator have different dispersion relations which can lead to an order of magnitude correction to 
leptogenesis in the weak washout regime and an order of 1 correction in the strong washout regime gSHlZ]. 

Keeping the above caveats in mind and applying finite temperature “cutting rules” (more discussion below), we 
obtain 


(A9) 



\M.^ (a —)■ zj) ^ — \M.^ (a —)■ zj) ^ 

rat 

'j' {T)c-aij{T) 

J2k,i 

\M° {a ^ kl)f + \M° {a ^ kl)f 

Cakl{T) 


where (a —>■ zj) is the amplitude for a —>■ ij at zero temperature and the sum over z'j' in the numerator is over 
intermediate states in the loop which go on shell as shown as the “cuts” in Fig. In Eq. (A9l the raij{T)'s are 
the thermal factors associated with the on-shell intermediate states, while the Caij(Tj’s are those associated with the 
final states. In the case of self-energy contributions, the factorized form as a product of thermal-dependent and zero 
temperature terms as in Eq. (A9 \ always holds (under the approximation that a is at rest with respect to the thermal 
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Figure 5. One-loop diagrams for decay a ^ ij. 


bath) while in the case of vertex diagrams, further approximations are required. One approximation we have made is 
to factorize out the temperature-dependent terms including the kinematic factors, and then to set the thermal masses 
in the rest of the terms to zero, which gives us expressions for these terms that coincide with the zero temperature 

results. In addition, we ignore the contributions from the cuts through a' and i', or a' and j' in the vertex diagrams, 

which as shown in Ref. |3S] in non-SUSY type-I leptogenesis can give corrections depending on the a — a' mass ratio, 
for example, at the level of 10% for ma'/ma = 1.1. Under these approximations, the temperature-dependent terms 
for both self-energy and vertex diagrams are the same and are given by 

Caij{T) = [1 -I- 77a (1 - 5ukj) + Vj^j)] A (1 -I- rjifP) (l -f VjfP) . (AlO) 

ra^](T) = [1 -I- 77a (1 - 6bi6bj) {rjiX^ + rjjXj)] A {l,Xi,Xj) (l -|- rjifP + VjfP) , (AH) 

with 6bi = 1 (0) if 7 is a boson (fermion) and 

A(l,a;,77) = \l{l+x-yf - Ax, 

E, = ^ (1 + Xi- Xj ), 


_ mi {Tf 
M 2 ’ 

Ej=Ma-E, = ^{l-x, + x,). (A12) 


For the statistical factors in raij{T), we applied the finite temperature cutting rules by considering causal (i.e., retarded 
or advanced) n-point functions as pointed out by Ref. |4H| which gives the dependence on the distribution functions 
^+TiifP+r]jfP in agreement with the results derived from nonequilibrium quantum field theory |48II501[62] in contrast 
to the results of Refs. [SUES] which obtained 1 -I- r]ifP -f Vj/P + ViVjfPfP when time-ordered n-point functions are 
considered instead. Notice that the imaginary time formalism also gives the statistical factors in agreement with the 
result of nonequilibrium quantum field theory gsHU]. 


Now we can apply the general results (AlOl and (All I to soft leptogenesis. For the decays N± —>■ and 

N± —>■ £Hu, the relevant thermal factors are obtained from Eqs. ( A10[) and (All) to be 


cf(T) = (l-xe- x^^^ A (1 - fP) (l - fpj , 

cb(T) = A (1,xj,xbJ (l + /p) (1 + fPJ , 

h(Y) = (^1-xg- x^^^ A (^l,xg,x^^^ (^1 - fP - jpJ , 

rsiT) = A {1,xj,xhP (l + /f + fP^) ■ (A13) 


The relevant thermal masses are m 


mj{Tf = 2mg{Tf = {^-gl + 

mnSTf = 2mj,pTf = Qgi + + ^A?^ T\ (A14) 


Next we prove a useful identity (in the context of the approximations made above) as follows 


rF{T)cB{T) - tb{T)cf{T) cx (l - fP - /£) (l + fp) (l + fPj - (l + /f + fP^) (1 - fP) (l - 


= - l) fPfp^e^«/^fpfP^ - - l) fpfP^e^-Z^fPjp^ 



= 0 . 


(A15) 
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In the second line above, we have made use of the following identity 

(1 + rjJP) (1 + rj.fp) = 


(A16) 


and the conservation of energy Ej:^ = Ei + E^ = Ej+ Efj^. Notice that this identity also holds if instead of using 
the factor 1 + rjif^ + rij/p in VatjiT), one uses 1 + + VjfP + ViVjfPfP as obtained in Refs. |3I1 IB5] . 

Finally, it can be shown that the CP asymmetries from gaugino contributions nanH] for the decays of N± to 
scalars and fermions are, respectively, given by egrF{T)cB{T) and —egrB(T)cF(T) where eg is some temperature- 
independent term. Using the identity (A15), these contributions sum up to zero illustrating the cancellations pointed 
out by Ref. pO] . 


Appendix B: Special cases of mixing CP asymmetries 

Here we will discuss two specific cases of mixing CP asymmetries. 

(a) Universal trilinear scenario: Aa = AYa- 

This is the scenario considered in Refs. [HITIIIS]. In this scenario, we are always in the regime (i) Ya ^ 
\Aa\/M, and from Eq. (18), we obtain [23] 


" cf{T) + cb{T) 


(Bl) 


with 


Im (A) ABTy 
M 4R2 + ■ 


(B2) 


(b) Simplified misaligned scenario: A„ = AY'^/(ZYa) ■ 


This is a specific scenario considered in Ref. m- In this scenario, we have from Eq. (|14|) 


G±(T) ^ ^ Y^^-^dcBin 


(B3) 


where d = ^ the minimum occurs at = 1/3 for all a. In Eq. (B3), we have dropped 

terms of 0{Y'^ms\jSY/M) except the second term which could dominate over the first when |Ap/M^ d~^. 
This condition can only be fulfilled if Pa deviates significantly from 1/3, i.e. a very hierarchical Pa- 

First let us consider the case |Ap/M^ <C d~^. From Eq. ([T7| we obtain 


4 < 


.Cf{T)-cb{T) 

- ■—FtAT- A) 

CF (T) -I- CB (T) 


1 |A| 


CB (T) 


9Pa 


1^1 _ J 

cf{T)+cb{tP^^ 


(B4) 


Regarding the first “thermal” term. Ref. m made a mistake in that the lepton flavors in the self-energy loop 
were not summed over resulting in expressions which were independent of Pa- Since this term coincides with 
Eq. (Bl I there is no enhancement for the simplified misaligned scenario case compared to the universal trilinear 


scenario case as claimed in Ref. 

In Ref. [TS], the second “nonthermal” term is ignored assuming it is always smaller than the first one. This is 
always true when the Pa’s are not hierarchical or only mildly hierarchical. The “nonthermal” term can dominate 
only in a hierarchical scenario when some of the Pa fulfill 


1 |A| 

P eg- _1 

“ ^ 3 M ■ 


(B5) 


Nonetheless the idea of considering a generic A term to enhance the efficiency was correct as shown in the current work. 
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When the condition above is fulfilled, we always have a mixed scenario where we have both thermal and 
nonthermal contributions - since Pa = 1, some of the Pq’s cannot fulfill Eq. (B5l. 


For the case |A| /M^ ^ d which can only happen when Pa is very hierarchical, we have from Eq. (171 




~ p 


M 2 _ cf{T)-cb{T) 

I .,2 , ——pjrT^—rB [T) 

\A\ d cb{T) 


1 

9 Pad 


CA Tb {T) 


(B 6 ) 


where 


_ Im (A) 

CA = 


dPEy 


M 


4P2 + ( 


lAl^ 

2M2 


The condition for the second “nonthermal” term in Eq. (B 61 to dominate is again Eq. (B51. 
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